Scanning Tunneling Spectroscopy on the novel superconductor CaCe 
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We present scanning tunneling microscopy and spectroscopy of the newly discovered supercon- 
ductor CaC6- The tunneling conductance spectra, measured between 3 K and 15 K, show a clear 
superconducting gap in the quasiparticle density of states. The gap function extracted from the 
spectra is in good agreement with the conventional BCS theory with A(0) = 1.6 ± 0.2 meV. The 
possibility of gap anisotropy and two-gap superconductivity is also discussed. In a magnetic field, 
direct imaging of the vortices allows to deduce a coherence length in the ab plane f a 6 ~ 33 nm. 



Interest in the superconducting properties of carbon- 
based compounds has been renewed by the discovery of 
superconductivity in CaC6 and YbC6 having a T c of 11.5 
K and 6.5 K respectively Q,Q|- Such values constitute 
an increase of nearly one order of magnitude of the T c re- 
ported for Graphite Intercalation Compounds (GIC) at 
ambient pressure, the latest being KTI1.5C4 with a T c 
of 2.7 K by Wachnik et a/.|3|. A controversy has arisen 
about the origin of superconductivity in these materials. 
On the basis of self-consistent electron-phonon calcula- 
tions, Calandra and Mauri |^j| suggested a conventional 
BCS phonon-mediated mechanism where the C out-of- 
plane and the Ca in-plane vibrations couple mainly to 
electrons of the Ca Fermi surface. A similar approach 
has been used by Mazin |5| with comparable results. On 
the other hand, Csanyi et aLp developed an excitonic 
model of superconductivity. Measuring the quasiparticle 
(QP) spectrum of these materials is thus of immediate 
interest. 

The preparation of high-quality samples by Emery 
et al. [7J has been an important step in the experi- 
mental study of CaCg : it has allowed a proper mea- 
surement of the bulk superconducting parameters J3| 
as well as the measurement of the magnetic penetra- 
tion depth, see Lamura et al.9}. Their results give evi- 
dence for a mainly s-wave superconductivity, and a BCS 
weak coupled superconductor with a gap A(0) = 1.79 
and 2A(0)/fcsr c =3.6 ± 0.2. However, in these mea- 
surements, the gap and the coherence length are ob- 
tained in an indirect manner. In this Letter, we report 
the first scanning tunneling microscopy and spectroscopy 
(STM/STS) measurements, as a function of temperature 
and magnetic field, giving directly the superconducting 
gap and the coherence length. Our results are in good 
agreement with a conventional BCS type superconductiv- 
ity. An upper bound for the gap anisotropy is estimated. 




FIG. 1: (Color online) (a) Topographic image of the cleaved 
surface of CaC6 (180 nmxl50 nm). The total height differ- 
ence is about 15 nm. (b) Temperature dependance of tunnel- 
ing conductance spectra. Experimental data (dots) are fitted 
using expressions (1) and (2) with V = 0.2 meV. The spectra 
are shifted by 0.5 for clarity, (c) Temperature dependance of 
the superconducting energy gap extracted from b). Solid line 
shows the BCS A(T) with T c =10 K and A = 1.6 meV. 



The samples were prepared by the immersion for ten 
days of a platelet of highly-oriented pyrolitic graphite 
in a molten lithium-calcium alloy at around 350 °C un- 
der argon atmosphere. This procedure, described in de- 
tail in |jj, leads to plate-like samples (2x2x0.3 mm in 
this study). They are polycristalline with c-axis of all 
the cristallites parallel to each other, whereas it is disor- 
dered in the ab plane. After the synthesis, the samples, 
still under argon atmosphere, were embedded into sil- 



ver epoxy between the STM sample holder and a cleav- 
age screw. This method allows to prevent air-sensitive 
compounds from damage prior to its introduction into 
the UHV chamber of the microscope. The samples were 
cleaved under UHV, and then rapidly cooled down to 4.2 
K, using a pure exchange He gas, at a pressure around 
1CP 3 mBar. Further cooling to 3 K was achieved by 
pumping on the helium bath of the cryostat. All the 
measurements were carried out using freshly cut Pt/Ir 
tips. In the chosen tip-sample configuration, tunneling is 
parallel to the c-axis. 

A typical STM topographic image is displayed figure 
Ula). One can find relatively flat regions, but the terraces 
are not as large as those usually encountered on precur- 
sor graphite samples. Also, it appears that the material 
does not cleave well. This illustrates that, despite the 
layered structure, CaCg is less anisotropic than most of 
the other GIC and has more pronounced 3-dimensional 
character, in agreement with Calandra et al. 4\. This 
is also confirmed by the coherence length reported by 
Emery et o/.|S| and by our measurement below. 

To obtain the conductance curves in the STS configu- 
ration, symmetrical bias voltage sweeps, typically in the 
range ±15 mV, were applied while acquiring the tunnel- 
ing current I{V). In a single measurement, 256 volt- 
age sweeps are applied to the junction (with feedback 
loop open) and the corresponding current-voltage spec- 
tra are averaged in order to eliminate the major part of 
the noise. The dynamical conductance curve dI/dV(V) 
is then the direct numerical derivative, without further 
data treatment. As is well known, at zero temperature, 
dI/dV(V) is directly proportional to the sample QP den- 
sity of states (DOS), Ns(Ep + eV). Taking into account 
thermal broadening leads to the expression : 



dI(V) 
dV 



dENs(E) 



where, 



N S {E) = N n (E F ) Re 



df{E - eV) 
dV 



E~iT 



y/(E-iT)*-A(ry< 



(1) 



(2) 



Here f(E) is the Fermi function, N n (Ep) is the normal 
DOS and T, in the BCS extended DOS, takes into ac- 
count the finite QP lifetime |lfl |. In these expressions, 
any Fermi surface or gap anisotropics are not taken into 
account. 

We performed tunneling spectroscopy in the tempera- 
ture range from 3.2K to 15K. The evolution of the con- 
ductance as a function of temperature is shown in figure 
IHb). The overall shape of the spectra agrees well with 
the extended BCS density of states (2). The set of curves 
can be fitted using a single isotropic gap A(T) and the 
temperature T as independant parameters, and T fixed 
at 0.2 meV in all the fits. The gap A(T) obtained from 
the fitting procedure is in good agreement with a con- 
ventional BCS self-consistent gap (figure Q]c), providing 



a zero-temperature gap value A(0) of 1.6 ± 0.2 meV and 
a T c of 10 ±1 K. This gives a BCS ratio 2A(0)/k B T c 
of 3.66, slightly higher than the expected value in the 
weak-coupling limit. T c and A(0) values are close to, but 
slightly smaller than, those previously reported: T c =11.5 
K by magnetization measurements and A(0)=1.79 meV 
obtained from magnetic penetration depth. The small 
discrepancy observed may originate from a little Ca de- 
pletion at the surface of the sample. 
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FIG. 2: (a) Effect of the parameter F on the conductance 
spectrum. Dots : experimental data at 3.2 K, dashed line: 
BCS fit with r=0, solid line: BCS fit with T= 0.2 meV. (b) 
Dots : experimental data, solid line: BCS fit with a cosine 
dependent anisotropic gap A=1.45 ± 0.55 meV. Inset: con- 
ductance spectra for different tunneling resistances (11 super- 
posed spectra from 0.1 Mfi to 10 MQ,). 

We now turn to a finer analysis of the tunneling spec- 
tra. In figure |2{a) we plot a typical dI/dV(V) exper- 
imental spectrum and two fits obtained by considering 
only a single isotropic BCS gap, at the temperature T = 
3.2 K of the junction. The first fit, with T set to zero, 
does not agree with the experimental data whereas the 
second fit with T = 0.2 meV is a perfect match. The need 
for this phenomenological smearing term, found often in 
the literature, raises the question of the real underly- 
ing physical phenomena. The T term was first suggested 
by Dynes 10], to take into account a finite QP lifetime 
(or imaginary part of the self-energy). However, in a 
STM experiment the broadening can be due to various 
effects : the finite QP lifetime corresponding to inelastic 
processes, the electronic noise of the STS, in particular a 
voltage jitter, and finally the anisotropy of the electronic 



structure, intrinsic to the material. For conventional su- 
perconductors, even strong coupling, the QP self-energy 
broadening is expected to be negligible. So, the use of a T 
parameter accounts for the smearing due to the intrinsic 
anistropy, and extrinsic experimental resolution, which 
must be distinguished. At least, a part of the T value is 
due to the resolution function of our experimental set-up. 

The possibility of gap anisotropy in this material is 
important to consider : the GIC's were believed to be an 
example of two-gap superconductivity 11]. This is due to 
the fact that the Fermi level is crossing both the graphene 
■k band, with a Fermi surface nearly cylindrical along the 
c* axis, and the intercalate s-band, with a Fermi surface 
almost spherical around the central T point of the Bril- 
louin zone. However, if the superconductor is anisotropic, 
then the tunneling geometry becomes important and one 
must express the precise tunneling of the sample QP, of 
wave vector k — (k\\,k±) and gap parameter Ag, to the 
tunneling tip, or the converse. 

Using the theory of Tersoff and Hamann [12J, N${E) 
in (1) must be replaced by the local tunneling DOS : 



N t{E,z) = ^J dS.^T^E^z) 



E 



E 2 - A 2 

k 



(3) 
where the integral is over the Fermi surface Sj? and the 
transmission factor is 13]: 



Tz(E,z) = \c u (E)\ 2 e- 2akZ 



(4) 



In (3), £fc is the normal spectrum, cr is the amplitude of 
the Bloch function at the surface of the sample, and a^ 

is the vacuum attenuation coefficient : a^ = J k 2 + ^p^, 
ip being the work function. Consequently, the main con- 
tribution to the tunneling conductance is due to states 
with small k\\ , or the surface Brillouin zone center (de- 
selection of the barrier). Neglecting the fc-dependence of 
T^(Ep, z) and taking a constant gap, gives back equation 
(2). The above equation (3) thus compactly represents 
the effects of Fermi surface, tunneling matrix clement 
and gap anisotropics. In the case of c-axis tunneling, the 
main contribution comes essentially from states of the 
spherical Fermi surface of the intercalant s-band. The 
inset of figure Hfb) presents the conductance spectra for 
different tunneling resistance, i.e. for different tip-sample 
distances. The spectra appear unmodified, providing ev- 
idence for the minor role of the e~ 2akZ term. In the case 
of a two-gap superconductor, one can observe additional 
signatures in the conductance. This can be understood 
by separating the integral (3) into different Fermi surface 
parts. Some of the contributions can be weak, if they are 
linked to Fermi surface points with larger k\\ . This situ- 
ation is observed in MgP>2 c-axis oriented thin films, the 
small gap A T dominates the spectrum but the QP peaks 
are also slightly affected. No double-gapped spectra, nor 
extra features, have been observed on CaCg. 



Since the precise value of (3) depends on the detailed 
band structure, in order to estimate a possible anisotropy 
we adopt a gap distribution of the form : 



A(0) = A + A x cos(20) 



(5) 



In this case, conductance (1) contains an additional in- 
tegration over the 6 parameter. The figure |2Jb) dis- 
plays a characteristic spectrum taken at 3.2 K and the fit 
with an anisotropic, cosine dependant gap with param- 
eters Ao=1.45 meV and Ai=0.55 meV. However, such 
a strong anisotropy is not supported by the fact that 
the spectra are unmodified with the different tip-sample 
distances (inset fig 2b) . Consequently, even if the occur- 
rence of a strong anisotropy in CaC6 or double-gap su- 
perconductivity cannot be completely excluded, it seems 
to be unlikely. Neverthless, a definitive answer would 
require tunneling along the ab direction. Such a mea- 
surement appears difficult, given the small thickness of 
the samples, and the air-sensitiveness of this compound, 
forbidding the ex-situ preparation of an inverted junc- 
tion, as was previously used on MgB 2 [lj] . It can also be 
reached with a very high energy resolution. This would 
require very low temperatures or the use of a supercon- 
ducting tip, as in the experiment of Rodrigo and Vieira 
|15j showing the multi-band superconductivity in NbSe2 . 

We now focus on scanning tunneling spectroscopy in 
the presence of a magnetic field. In a ideal type II su- 
perconductor, an applied magnetic field B a penetrates 
the sample in the form of vortices, each carrying one 
flux quantum <I>o. Due to energy considerations, vortices 
arrange in a periodic lattice, usually triangular with a 
spacing d = (2 V3 Qa/Ba) 1 / 2 . Each vortex is surrounded 
by screening currents which decay over the magnetic pen- 
etration length A and has a core extending over the co- 
herence length £, where superconductivity is suppressed. 
Finally, when the thermodynamic critical field H c i is 
reached, the vortices overlap and the sample return to 
normal state, with a metal-like QP DOS. 

In order to determine £, we have attempted to image 
the vortices by scanning tunnelling spectroscopy under a 
magnetic field applied perpendicular to the sample, thus 
to its c-axis. This configuration enables to determine 
only the in-plane coherence length £„&, for the reasons 
of selective tunneling discussed earlier. As also men- 
tioned, the samples do not cleave very well, and it has 
been difficult to obtained surfaces suitable for scanning 
spectroscopy. Several attempts were necessary in order to 
image the vortices, which appear clearly only in limited 
regions. The inset of figure [3] a) displays the normalized 
zero bias conductance map, obtained for a magnetic field 
S a =0.05T where bright colors correspond to vortex cores 
and dark colors to superconducting areas. For this value 
of B a , the inter- vortex distance is of the same order than 
the magnetic penetration depth A |2(. In this case, the 
screening currents are expected to affect the QP DOS in 
the superconducting region between the vortices. This 
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FIG. 3: (Color online) (a) Normalized zero bias conductance 
versus position from the vortex center. Inset: Zero-bias con- 
ductance map showing three vortices (170 nnx230 nm). (b) 
Evolution of the dl /dV(V) spectra on a grey scale as a func- 
tion of the distance to the vortex core, (c) Normalized con- 
ductance spectra with increasing magnetic field at T=5.5K. 



can be seen in the figure EJb) which displays the evo- 
lution of the dI/dV(V) spectra on a grey scale as one 
goes from the inter-vortex space to the core. Outside the 
core, the zero bias conductance is significant and the QP 
peaks are affected. A normalized zero bias vortex profile 
is show figure|^a). Following reference [lj(, the zero bias 
profile can be fitted by the formula below, derived from 
the Ginzburg-Landau expression for the superconducting 
order parameter : 

a(r, 0) = ao + (1 - cto) X (1 - tanh(r / (V2(,)) (6) 

where &q is the normalized zero bias conductance away 
from a vortex core and r the distance to the vortex center. 
The fit yields a coherence length of £ a fc=46 nm at T=5.5K 
and an extrapolated £ Q {,(0)=33 nm at zero temperature. 
This value is in a good agreement with £ o t,(0)=35 nm 
reported from bulk magnetic measurements [8j- 

In the case of a clean superconductor, Andreev bound 
states [la| due to the confinement of QPs inside the vor- 
tex core are expected to appear as a zero bias conduc- 
tance peak in the tunneling conductance |l7j . However, 
in the dirty limit, the scattering of these states results in 
a smearing of the peaks to bumps and, eventually, to a 
flat metal-like DOS. This has been previously observed by 
STS measurement in Nbi_ x Ta a ,Se2 at 1.3 K |ljj. Even 
if some of our spectra display some weak sub-gap fea- 
tures close to the zero bias, similar to those of reference 
|l8| , there is no clear evidence for the presence of bound 



states. This is consistent with the fact that these sam- 
ples are in the dirty limit, in accordance with ref |9|]. The 
figure Bib) displays the evolution of conductance spectra 
at 5.5K with further increase of magnetic field. The main 
effect is to smear the QP peaks and to fill the gap pro- 
gressively with QP states. For a field B a ~ 0.6 T, the 
conductance displays a normal DOS indicating that the 
sample has returned to the normal state. This value is 
larger than the one extracted both from magnetometry 
measurement and from the value of £ deduced above. In 
this "zero-field cooled" experiment, this is likely due to 
flux trapped in some poorer quality area of the sample 
which modifies locally the magnetic field distribution. 

In conclusion, we have performed tunneling spcctro- 
copy measurements on the superconducting Ca interca- 
lated graphite, CaCg, at various temperatures. The spec- 
tra and their temperature dependence are consistent with 
a single isotropic BCS gap of 1.6 meV, with no clear indi- 
cation for additional contribution. Nevertheless, further 
measurements in different tip-sample configurations (a- 
axis for instance) will be necessary to give a definitive 
conclusion concerning the possibility of two-gap super- 
conductivity or strong gap anisotropy. In the presence of 
a low magnetic field, the coherence length £ a ;, = 33 nm 
is extracted directly from the real space vortex imaging. 

This work was supported by project GPB Materiaux aux 
proprietes remarquables. 



[3 

[4; 

[5] 
[6] 

[7] 

[s; 

[9 
[10 

[11 
[12 

[13 
[14 
[15 
[16 

[17 
[18 
[19 



T. E. Weller et al., Nature Physics 1, 39-41 (2005). 

Of historical interest, CaC6 was first prepared in 1980 (D. 

Guerard, M. Chaabouni. P. Lagrange. M. El Makrini and 

A. Herold, Carbon 18 257-264 (1980)). Just like MgB 2 , 

CaC6 is "an old material but a new superconductor" . 

R. A. Wachnik et al., Solid state commun. 43, 5-8 (1982). 

M. Calandra, F. Mauri Phys. Rev. Lett. 95, 237002 

(2005). 

I. I. Mazin Phys. Rev. Lett. 95, 227001 (2005). 

G. Csanyi et al., Nature Physics 1, 42-45 (2005). 

N. Emery, C. Herold, Ph. Lagrange, J. Solid State Chem. 

178, 2947-2952 (2005). 

N. Emery et al., Phys. Rev. Lett. 95, 087003 (2005). 

G. Lamura et al., Phys. Rev. Lett. 96, 107008 (2006). 

R. C. Dynes, V. Narayanamurti and J. P. Garno Phys. 

Rev. Lett. 41, 1509 (1978). 

R. Al-Jishi, Phys. Rev. B 28, 112-116 (1983). 

J. Tersoff D. R. Hammam Phys. Rev. Lett. 50, 1998 

(1983). 

Higher Bloch terms have been neglected. 

F. Giubileo et al., Phys. Rev. Lett. 62, 177008 (2001). 

J. G. Rodrigo, S. Vieira, Physica C 404, 306, (2004). 

C. Caroli, P. G. De Gennes, J. Matricon Phys. Lett. 9, 

307 (1964). 

H. F. Hess et al. Phys. Rev. Lett. 62, 214 (1989). 

Ch. Renner et al., Phys. Rev Lett. 67, 1650 (1991) 

M. R. Eskildsen et al. Phys. Rev. Lett. 89, 187003 (2002) 



